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In homogeneous and isotropic loop quantum cosmology, gravity can behave repulsively at Planck- 
ian energy densities leading to the replacement of the big bang singularity with a big bounce. Yet in 
any bouncing scenario it is important to include non-linear effects from anisotropies which typically 
grow during the collapsing phase. We investigate the dynamics of a Bianchi I anisotropic model 
within the framework of loop quantum cosmology. Using effective semi-classical equations of motion 
to study the dynamics, we show that the big bounce is still predicted with only differences in detail 
arising from the inclusion of anisotropies. We show that the anisotropic shear term grows during 
the collapsing phase, but remains finite through the bounce. Immediately following the bounce, the 
anisotropies decay and with the inclusion of matter with equation of state w < +1, the universe 
isotropizes in the expanding phase. 
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I. INTRODUCTION 

Recently, the investigation of loop quantum cosmol- 
ogy (For a review see jli) in homogeneous and isotropic 
universes has indicated that the classical big-baiig singu- 
larity can be replaced with a big-bounce [3, H, 0, S @] ■ 
In these scenarios gravity can be interpreted as becom- 
ing repulsive in the Planckian high energy regime, im- 
plying that our current expanding universe would have 
been preceded by a contracting phase. The presence of 
a big-bounce has been shown to be a rather generic and 
genuine feature of the quantum gravitational effects of 
loop quantum cosmology and does not require any exotic 
matter which violates energy conditions [7|. An exciting 
fact of bouncing cosmological models is that the scales 
measured in the cosmic microwave background (CMB) 
can be in causal contact if the current expanding phase 
is preceded by a contracting one, thus opening the possi- 
bility for a replacement of the standard inflationary sce- 
nario. 

In order to more fully develop this scenario, one must 
go beyond the assumption of homogeneity and isotropy. 
It is both the inhomogeneities and anisotropies that are 
expected to grow in a collapsing phase and thus a proper 
accounting of these fluctuations is required. One par- 
ticular question that immediately arises is whether the 
presence of the bounce is stable under the inclusion of in- 
homogeneities and anisotropies. A proper description of 
the inhomogeneous perturbations and anisotropies would 
then provide an answer to the question of whether a suit- 
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able alternative to inflation can be constructed, and/or 
what possible cosmological signatures may result. 

In this paper we do not seek to answer all these ques- 
tions as work on including inhomogeneities in loop quan- 
tum cosmology is in its infancy (initial progress can be 
found in [1, 0, [l^, [HI). Instead, we will focus on the 
behavior of anisotropies, studying the dynamics of the 
anisotropic Bianchi I model in the framework of loop 
quantum cosmology. In the classical Bianchi I universe 
sourced with matter with zero anisotropic stress, the 
anisotropic shear term behaves as an effective matter 
component with energy density that scales as in the 
Friedmann equation with a being the mean scale fac- 
tor. Thus for matter with equation of state w < +1, the 
anisotropies will dominate the collapsing phase as the sin- 
gularity is approached. If the current expanding phase 
of the universe was preceded by a collapsing phase, the 
inclusion of anisotropies can be expected to play a signif- 
icant result near the bounce. 

The loop quantization of the anisotropic Bianchi I 
model was initially studied in \12\ and more recently 
in [Tsj . In this paper, we will study the dynamics of 
the model at the level of effective classical equations of 
motion that incorporate quantum effects arising from 
the loop quantum Einstein equations. At the level of 
the effective equations we study, we show that the big- 
bounce is indeed robust under the inclusion of non-linear 
anisotropies. We will show that the anisotropic shear 
remains finite through the bounce and that if matter 
with equation of state < +1 is included, the universe 
isotropizes in the expanding phase. The results represent 
evidence that the bouncing scenario of loop quantum cos- 
mology is robust when the assumptions of homogeneity 
and isotropy are relaxed and gives hope that the same can 
be said when inhomogeneous perturbations are properly 
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included in the theory. 



II. CLASSICAL DYNAMICS 

We start with the classical setup for the Bianchi I 
model. Loop quantum cosmology is based on a Hamil- 
tonian formulation and thus we will define the Hamilto- 
nian for the gravitational plus matter degrees of freedom. 
Hamilton's equations of motion then are equivalent to 
Einstein's equations for the model considered. Since the 
starting point of the quantization is the Hamiltonian, it 
is there that we will incorporate the quantum effects in 
a modified effective Hamiltonian to be introduced in the 
next section. 

Starting with the dynamical variables that comprise 
the classical phase space, we have three triad variables pi 
and three connection variables c/ with / = 1,2,3. The 
classical metric given in terms of the directional scale 
factors a/ is 



'N'^ dt^ + a\ dx^ + aj dy'^ + al dz^ 



(1) 



with N representing the lapse which is a freely specified 
function representing the freedom to redefine the time 
variable, and the coordinates x, y, z are all valued on the 
entire real line (we are not considering a compactified 
model). The triad variables are directly related to the 
scale factors as 



|pi|=a2a3, |p2| = aias, IpsI = aia2. 



(2) 



Thus the triad variables encode information about the 
spatial geometry. The connection c/ will encode infor- 
mation about the curvature (essentially time derivatives 
of the scale factors) as will be evident once Hamilton's 
equations are solved. 

In the Hamiltonian framework, the classical equations 
of motion are derived from a Hamiltonian obtained by 
inserting the homogeneous phase space variables into the 
Hamiltonian of general relativity making it a functional 
of p/ and c/ in this case. However, a non-triviality arises 
from the spatial integrations in the Hamiltonian. Be- 
cause of homogeneity, these integrations diverge since 
we are considering the non-compact Bianchi I model. 
To overcome this, the spatial integrations can be re- 
stricted to a finite sized fiducial cell with fiducial volume 
Vo = J d^x. We can use the fiducial cell to define untilded 
variables pi , c/ as 



P/ = Vg ' PI 



Cl 



T/l/3~ 

Vn' CI. 



(3) 



In the classical theory, one can freely rescale the coordi- 
nates x,y,z while leaving the physics invariant. Under 
this rescaling, one can show that the untilded variables 
Pi,ci are invariant. Since the quantization is based on 
the untilded variables, the quantum theory is manifestly 
invariant under this coordinate rescaling. This is entirely 
analogous to the non-compact isotropic fc = [l3| and 
A; = — 1 models where the same procedure is used. 



With the understanding that the spatial integrations 
in the Hamiltonian are restricted to the fiducial cell, the 
total Hamiltonian of the model is given by 



n 



^giav H~ '^matter 



-N 



+'Hm ■ 



(C2P2C3P3 + C1P1C3P3-I-C1P1C2P2) 



(4) 



with TLm being the matter Hamiltonian. Here we have 
K = SttG, and 7 is known as the Barbero-Immirzi param- 
eter and represents a quantum ambiguity of loop quan- 
tum gravity which is a non-negative real valued param- 
eter. Two sets of equations of motion derived from the 
Hamiltonian then govern the dynamics. First, the Hamil- 
tonian of gravity and matter is of the constrained type 
whereby it vanishes identically for solutions to Einstein's 
equations. Thus an equation of motion is given by 



n = Q. 



(5) 



Second, Hamilton's equations give the time evolution of 
any phase space variable through the Poisson bracket: 



P'/ = {p/,^}, c/ = {c/,7Y}. 



(6) 



The Poisson structure of the gravitational variables leads 
to the only non- vanishing Poisson brackets 



(7) 



With this, Hamilton's equations for c/ and pi are given 
by 



an 



(8) 



To proceed further with the equations of motion, we 
must specify the matter Hamiltonian. We take it to be 
of form 



Ti-M = N^pip2P3 PM 



(9) 



with PM being the matter energy density. In this paper, 
we assume that the matter has zero anisotropic stress 
which implies that pM couples to pi in the form 



Pm{P1,P2,P3) = PMiPlP2P3) ■ 



(10) 



This assumption is true for scalar fields and perfect fiuids 
which we will concentrate on in this paper. Additionally, 
to derive the equations of motion we must specify the 
form of the lapse. For simplicity of the resulting equa- 
tions, let us choose a form given by 



(11) 



which means we are using a different time t' variable than 
the usual cosmic time t given by 



dt' = {piP2P3)'^'^dt. 



(12) 
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Note that the choice of lapse is arbitrary and none of the 
physical results depend on the choice. 

We now derive the classical equations of motion. The 
Hamiltonian with our choice of lapse is 



n 



1 



(C2C3P2P3 + C1C3P1P3 + C1C2P1P2) 



ViViVz Pm- 



(13) 



The first set of Hamilton's equations in ([8]) for the time 
evolution of the triad then yields for instance 



dpi Pi , , . 



dt' 
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(14) 



Using the relations between the triad components and 
scale factors ([2]) and ([3]) , we can solve these equations for 
the connection coefficients c/ to get 



c/ = 7K) (010203) —=-iVq' — 



dt 



(15) 



These relations therefore provide us with the interpreta- 
tion of the connection components as containing infor- 
mation about the curvature, which in the the Bianchi I 
model is entirely encoded as the extrinsic curvature here 
proportional to daj/dt. 

The next set of Hamilton's equations in ([5]) for the time 
evolution of the connection yields for instance 



dci 
IF 



(C2P2 + C3P3) 

7 



+K7P2P3 Pm +P1 



dpM 
dpi 



Combining (fH|) and (fTO|) gives a key relation 
d , . ( ^ dpM 



(16) 



(17) 



It is here we use our assumption that the matter has 
zero anisotropic stress (jlOp . This assumption implies 
PidpAi/dpi — pjdpM/dpj and thus p7)) yields 



— {pid -pjcj) = 0, 

which can be integrated to give 

Pici - pjcj = aij 



(18) 



(19) 



with aij being a constant anti-symmetric 3x3 matrix 
satisfying by construction ai2 + 023 -|- a^i — and the 
factors of 7, Vq are chosen for convenience. Written in 
terms of the scale factors this implies 



Hi - H.J = 



CtlJ 

010203 



(20) 



for the directional Hubble rates (in terms of cosmic time 
dt = ^9010203 dt' ) 



Hj 



21 
ai 



(21) 



Using pO]) and the vanishing of the Hamiltonian, we 
can write down a generalized Fricdmann equation. The 
vanishing of the Hamiltonian (|13p gives in terms of the 
Hubble rates 

H,H2 + H,H3 + H^H-s = npM • (22) 

Let us next define the mean scale factor a as 

a = (010203)1/3, (23) 

from which the following relation holds 

o _ 1 
a 3 

This relation implies 

2 



(24) 
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1 

~ -^\HiH2 + H1H3 + ^2^3) 
[iHi-H2f + [Hi-H-if + {H2-H:if]. (25) 



Finally, using (PU)) and (P^ . the generalized Friedmann 
equation becomes 



■ \ 2 



at 6 0° 



(26) 



where S is given in terms of the constants of motion a/j 
as 



= Y8 '"12 + "23 + "31 J 



(27) 



The anisotropic shear cr^^cr'^'' is related to S through 



1 



3 
6E2 



{Hi-H2f + {H2-H3f 



(28) 



The equations of motion for a(t) can now be determined 
from the generalized Friedmann equation (I26p . The time 
evolution for the directional scale factors ai{t) can then 
be determined from a combination of ((20l) and ((24| , once 
a{t) is known. 

Let us discuss the interpretation of the generalized 
Friedmann equation . We first recognize the isotropic 
matter term k/9m/3 on the right hand side with the 
anisotropic shear term behaving as a stiff fluid with equa- 
tion of state w = +1. The isotropic limit is achieved 
when aij = 0, E = which from the relations (|^D)) im- 
plies that the directional Hubble rates Hj are identical 
indicating isotropic expansion. Because the shear terms 
scales as o^^, for typical forms of matter the early uni- 
verse can be anisotropy dominated, with matter then 
dominating the later stages. In particular, if w < -1-1, 
the later stages of expansion tend to a more isotropic 
state, whereby the Hubble rates become identical in all 
directions. Oppositely, as the singularity is approached, 
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the universe behaves hke a vacuum universe as the shear 
term dominates the matter term and Kasner hke behav- 
ior occurs. This behavior is desribed in more detail in 
appendix \^ 

As an example, let us consider a dust filled (w = 0) 
Bianchi I universe with matter density pM — A/a^. The 
generalized Friedmann equation can be solved analyti- 
cally giving 



a{t) 



3kA 



1/3 



(29) 



which gives the standard dust filled isotropic behavior 
at late times a cx t^/"^ while at early times is anisotropy 
dominated a cx t^^^ and is singular at i = 0. The time 
evolution of the individual scale factors can then be de- 
termined from (|20|) and ([24)) . which gives for instance 



ai{t) = 



3S 



^12+"31 



i{t). 



(30) 



Let us consider the situation where initially ai is con- 
tracting with the other two directions expanding. In 
this case we have ai2 < and a^i > 0. It is easy 
to see from ([30]). that at late times ai{t) behaves as 
ai(t) « a(t) cx i^/^ and thus is expanding in accordance 
with the isotropic behavior. At early times, however, it 
is not too difficult to show that pO|) implies contraction 
(di < 0) and thus the scale factor bounces at some finite 
time and eventually the universe isotropizes to a matter 
dominated phase. Therefore with matter with equation 
of state w < +1, the universe can initially be in a Kasner 
like epoch with one direction contracting and the other 
two expanding; eventually, however, such a universe will 
isotropize giving expansion in all three directions. 

Thus, at the classical level the anisotropics are ex- 
pected to play an important role as the singularity is 
approached. In the isotropic setting of loop quantum cos- 
mology, the bouncing scenario has been well described in 
terms of an effective Friedmann equation of the form [3] 



rPAf(l 



Pm . 

Pcrit ' 



(31) 



with pcrit ~ ■82ppi being a critical energy density of 
Planckian order. When the matter density reaches the 
critical density, H — indicating a bounce. Because the 
bounce occurs in the high energy regime near the classi- 
cal singularity, it is important to include the effects of the 
anisotropics. In the next section we will answer the ques- 
tion as to what role the anisotropics play in the bouncing 



scenario. 



III. EFFECTIVE LOOP QUANTUM DYNAMICS 

As we have stated, the quantum modifications due to 
loop quantum cosmology that we study exhibit them- 
selves in the form of a modified effective Hamiltonian. 



Since we study effective classical equations of motion, 
we are by definition ignoring certain quantum degrees 
of freedom. Properly speaking, in the quantum the- 
ory, dynamics is understood through expectation values 
of observables calculated from semi-classical wave pack- 
ets. Thus additional effects which we do not study can 
arise from features of the semi-classical state such as the 
spread and so forth. The effective quantum modifications 
we study to first order are insensitive to the features of 
the semi-classical wavefunction and are expected to be 
valid provided the wavepacket remains sharply peaked 
though the evolution. That this is a good approxima- 
tion has been verified in the isotropic models of LQC 
sourced with a massless scalar fieldP, d, @, where the 
quantum dynamics have been extensively developed and 
understood. Note that the effective equations we consider 
are somewhat heuristically motivated and more system- 
atic approaches to deriving effective equations have been 
considered [15,] , but ultimately proper justification of any 
effective scheme requires the study of the quantum dy- 
namics. 

The main quantum effect arises from the fact that 
in LQC, the connection variables c/ do not have di- 
rect quantum analogues and are replaced by holonomies 
(roughly exponentials of the connection). This manifests 
itself in the effective Hamiltonian by replacing the clas- 
sical c/ terms with sine functions 



c/ 



(32) 



where flj are real valued functions of the triad coefficients 
Pi which are a measure of the discreteness in the quantum 
theory. Classical behavior is expected in the limit when 
/i/C/ <^ 1 whence sin{fljci)/ p,i « c/. 

In the original construction of LQC, the /2/ were taken 
to be constants (referred to as ^o)llJ|i however it has 
been shown that in the isotropic case this can lead to 
the wrong semi-classical limit [J]. In the isotropic case it 
has been argued that fl should not be a constant, but 
should scale as /l cx l/^/p which has been shown to have 
a nice semi-classical limit [4]. Extending this scheme to 
the Bianchi I model is slightly more ambiguous as several 
possibilities exist. We choose in this paper to focus on 
the scheme proposed in |13|| where /2/ are given by 



= 1 — 



(33) 



where A — ^{4T:^£p^) is the area gap in the full theory 
of LQC and ipi = \J Gh is the Planck length. An alter- 
native scheme would be to have fl i oc 1 /aj which we will 
briefly discuss in Appendix [C] However the main physi- 
cal results we will describe do not depend sensitively on 
either of the scheme; only the quantitative results would 
change. Better input from the full theory might pro- 
vide justification for either scheme. One advantage of 
the scheme ([55]) is that it is more amenable to the study 
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of dynamics with semi-classical states in the quantum 
theory, as the resulting difference equation is much sim- 
pler. This would allow for the test of the validity of the 
results presented here, by examining semi-classical state 
behavior in the quantum theory. 

Additional modifications due to loop quantum cosmol- 
ogy have been studied extensively in the isotropic setting 
and pertain to operator eigenvalues of inverse pi factors 
that would appear in the matter part of the effective 
Hamiltonian. In the isotropic setting this amounts to re- 
placing factors of p~^^^ in the matter Hamiltonian by an 
operator eigenvalue function dj (p) which is bounded and 
vanishes at p = corresponding to the classical singu- 
larity.^ Additionally an ambiguity parameter j appears 
such that larger values amplify the dj effects. In this 
paper we will for simplicity ignore these effects in ac- 
cordance with arguments that j should take its smallest 
value [H, [3] J as well as questions as to the ambiguous 
nature of the critical scale at which the corrections are ap- 
preciable. Further discussion can be found in [3. Isl. M. . 
Our results presented here will remain valid as long as the 
critical scale remains below the scale at which the bounce 
occurs. In lieu of this, in our analysis the matter energy 
density pM appearing in the effective Hamiltonian p4|) 
will not contain factors of dj and will assume the classical 
form. 

With these caveats in mind, the effective Hamiltonian 
with lapse N — y^piPiPs is given by 



Ti-cs = 7 



1 rsin(^2C2)sin(^3C3) 



■P2P3 



l_ /i2^3 

-l-cyclic terms \ + P1P2P3 Pm (34) 



and it becomes easy to see that in the limit of small p-icj, 
the classical Hamiltonian ((T3]) is recovered. 

Hamilton's equations proceed in the same fashion as 
the classical setup. The equations for dpi/dt' and dci / dt' 
give for instance 

dpi piCos(/IiCi) r sin(^2C2) , sin(^3C3)\ 

-777 = 1^2 : +P3 = \ 35) 

at' 1 i P2 P3 ) 



and 



dci 1 /3sin(/iiCi) ciCOs(/iiCi) 



dt' 7 V 2^ 



sm(/i2C2) sm(^3C3) 
X <,P2 : hp3 : 

P2 P3 

dpM ' 
dpi 



+ K7 P2P3 PM + Pi 



(36) 



Hamilton's equations are now more complicated since the 
discreteness parameters /2/ depend on pi, but let us de- 
fine Gi{t') as 



Giit') -.^pi 



sin(/i/c/) 



PI 



(37) 



which can be shown to satisfy 



^ = K7C0s(^/C/)piP2P3 i^P^'^ ^ ■ (^^) 

The vanishing of the Hamiltonian gives 

ni^PiP2P3 Pm = Q1Q2 + Q1Q3 + Q2Q3- (39) 

Note that the classical limit is attained by /I/C/ 0, 
where we have sin(/2/C/)//i/ cj, cos{p,jCj) 1 and 
therefore ([55]) . ([55)1 reduce to their classical coun- 

terparts (Ull), (dll), pT]). 

At this stage, the equations are too complicated to 
solve analytically as was possible classically. To get a 
handle for the evolution, we can consider as an example, 
the case of a massless scalar field with equation of state 
w = whereby the equations simplify considerably. 
With a massless scalar field, we have 



Pm 



2piP2P3 



(40) 



where the momentum P? is a constant of motion. There- 



fore, the term pM + Pi 



dpM 
dpi 



vanishes identically. 



From this wc find that Qj are all constants in time. Equa- 
tion (|37p then implies that the triad components arc all 
bounded as 



PI > [\gi\VA 



2/3 



(41) 



implying that the classical singularity is never ap- 
proached. One can show from the equations of motion, 

that Pi > when pi = (|C//|\/A)^^'^ implying that the 
individual triad components pi bounce and hence the 
whole universe must bounce. With this bound one can 
show that no curvature invariants blow up, and the dy- 
namics is non-singular.^ Note that this analysis holds for 
the vacuum case if P^ = 0. Therefore a bounce would 
occur generically for the vacuum case also. 

Thus, at least with a massless scalar field, the bounce 
is robust under the inclusion of anisotropics. Next we 
can monitor the behavior of the anisotropic shear term 
through the bounce for this example. If we look at the 



^ Similar inverse triad effects appear in the gravitational Hamil- 
tonian leading to a function Sj{p) [IflL These corrections have 
been typically been ignored, but see |17j for work on some cos- 
mological implications. 



^ There is a special case when one of the Qj vanishes whence pj 
is not constrained above a minimum value. However, this is the 
special case when pi is a constant which can be seen from 1351 1. 
Therefore the specific triad component does not become zero and 
the evolution remains non-singular. 
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shear parameter E we showed that classicany it remains 
constant throughout the evolution. With the loop quan- 
tum modifications we can monitor its behavior by using 
its definition 

■■=^(^{Hi-H2)^ + iH2-H3f + iH3-Hi)^y (42) 

which in the classical case was a constant because of equa- 
tion pTjl . In the LQC case, let us assume we start with a 
nearly classical contracting universe with each /1/C/ <C 1. 
In this limit from ([37]) we find 

Gi ~ Pici. (43) 

Since Gi are constant for the massless scalar field, we can 
identify them with the classical constants q/j from ([TO]) 
as Qi — Qj = iV^aij. The shear factor S is then given 
initially as 

S^(pre bounce) « —^-rj- { {Gi - G2Y + cychc terms ) . 

1V>^Vq \ ) 

(44) 

After the bounce, the constancy of Gi and equation (137)) 
imply that /2/C/ become small as p/ grow and hence clas- 
sical behavior is recovered at late times. We can apply 
the same argument to conclude that the late time behav- 
ior of S approaches the pre-bounce value in terms of the 
same constants Gv- 

S^(post bounce) = S^(pre bounce) (45) 

and hence the shear factor E is conserved before and 
after the bounce. Using equations ([35)1 and ([37)) . it is not 
too difficult to show that the shear term remains finite 
through the bounce indicating that the anisotropics do 
not blow up. 

The similar conclusions can be obtained for the generic 
cases with the inclusion of arbitrary matter with w < -1-1. 
The details of the effective loop quantum dynamics with 
generic perfect fluids are investigated in Appendix |B] 
The detailed analysis shows that individual -pi in differ- 
ent directions can bounce at slightly different moments. 
Moreover, there is competition between the matter en- 
ergy density pu and the directional density gj (defined in 
(|B10[) . which is associated with the classical anisotropic 
shear) to be dominant when the quantum corrections 
start to become significant. As a result, depending on 
how energetic the matter content is (compared to the de- 
gree of anisotropics), the bounce can take place either 
in the "Kasner phase" or in the "isotropized phase" (or 
in the "transition phase" in between). If quantum cor- 
rections take effect in the isotropized phase, the bounce 
happens around the moment when pM approaches .82ppi, 
giving similar results as in the isotropic model. On the 
other hand, if the bounce occurs in the Kasner phase, the 
individual triad components bounce when gj approach 
.86ppi. 

In the case that the big bounce occurs in the 
isotropized phase, the Kasner fashion of the classical so- 
lution is smeared by the quantum effect and thus the 



information of anisotropics is blurred. Therefore, on 
the other side of the bounce, the classical anisotropic 
shear is changed. This explains why E^(post bounce) 7^ 
E^(pre bounce) in general. Only in the case when the 
big bounce takes place in the Kasner phase do we have 
S^(post bounce) ~ E^(pre bounce)."^ 

In the next section we explicitly show the results men- 
tioned here by numerically solving the equations of mo- 
tion. 



IV. NUMERICAL RESULTS 

Owing to the complexity of the effective equations of 
motion, we can not go further in the analytical analysis 
for more general forms of matter. In this section we will 
present some numerical simulations of the equations of 
motion by including other forms of matter. In particular 
we will focus on whether a bounce occurs for an initially 
contracting universe and how the shear term E behaves 
through the bounce. Additionally we will show that the 
universe isotropizes after the bounce with the inclusion 
of w < -1-1 matter. 

The differential equations for the time evolution of 
Pi{t') and cj{t') are given in ([35]) and ([35)1 . Numerically 
we solve the equations given a lapse N equal to one and 
plot as a function of cosmic time t. The initial condi- 
tions are chosen consisting of a collapsing semi-classical 
universe such that d < and flicj <C 1. The shear term 
E is monitored using the classical formula (pS)) . 

The first example we consider is the vacuum case. The 
analysis can be understood analytically as a special case 
of the massless scalar field studied in the previous section. 
The same analysis holds when = which gives the 
vacuum case. Therefore, for a vacuum case a bounce is 
generic and the shear term E^ is conserved before and 
after the bounce. 

The mean scale factor a{t), directional scale factors 
ai{t), and shear term E^(<) are plotted in Figures [1] and 
[2] respectively, for a representative numerical simulation 
consisting of an initially contracting universe which in the 
vacuum Kasner case requires two contracting directions 
and one expanding. After the bounce, the universe con- 
sists of a Kasner like expanding phase with growth in two 
directions and contraction in the other. Throughout the 
evolution, the shear term remains finite with the behav- 
ior during the bounce showing no obvious pattern (from 
numerical simulations, the exact behavior can depend on 
initial conditions and choice of matter). However, it is 
clear that the post-bounce value tends to the pre-bounce 



In the cases of vacuum and with only massless scalar field, the 
classical solution remains in Kasner phase throughout the evo- 
lution (the universe is not isotropized by matter). Thus, quan- 
tum corrections always take effect in the Kasner phase and this 
is why we have S^(post bounce) = S^(pre bounce) exactly for 
these two special cases. 
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value as shown by the analytical arguments of the previ- 
ous section. 

The inclusion of a massless scalar field does not change 
the results of the vacuum case significantly. With matter, 
the initial contracting phase can be either one where all 
directions contract, or a Kasner like contraction as in 
the vacuum case. In Figures [3] and [H the scale factors 
and shear term are plotted for the case where all three 
directions initially are contracting. Again, the shear term 
is finite through the evolution and the post-bounce value 
approaches the pre-bounce value, as expected from the 
analytical analysis. 

As an example of the inclusion of other forms of mat- 
ter, we consider a radiation (w = 1/3) field with energy 
density 

PM = ^ (46) 

with A some constant. The numerical behavior of the 
scale factor appears in Figure [5l where again the bounce 
occurs. Various initial conditions were examined and in 
all cases the bounce occurred providing evidence that it 
is a general feature of the effective Hamiltonian we study 
in this paper. In contrast to the massless scalar field 
case, the shear term is not conserved after the bounce. 
This is shown in Figure [51 and for the particular initial 
conditions chosen, actually decreases in the post-bounce 
regime. It must be stated that the decrease is not a 
generic feature and is sensitive to the initial conditions. 
However, as in the previous cases, the shear term is finite 
through the entire evolution. 

Since the equation of state for the radiation field w — 
1/3 < +\, the post-bounce regime leads to an isotropiza- 
tion of the universe. This is borne out in the numeri- 
cal simulations in Figure [71 There are plotted ratios of 
the directional Hubble ratios. The ratios approach unity 
in the post-bounce epoch indicating identical expansion 
rates in the three directions. 



V. DISCUSSION 

Let us restate the main results presented. We have con- 
sider the anisotropic Bianchi I model with loop quantum 
corrections to the classical equations of motion. We have 
shown that a bounce occurs under rather generic con- 
ditions in a collapsing universe. The anisotropic shear 
term S^/a^, which classically grows during collapse and 
blows up at the singularity, remains finite through the 
bounce. After the bounce, the universe behaves more 
and more classically and can isotropize at later times. 
This is thus evidence that the bouncing scenario of the 
isotropic models of loop quantum cosmology is robust 
when the symmetries of the isotropic model are relaxed. 

We must reiterate the caveats that has gone into this 
analysis. First, we have used entirely the effective clas- 
sical equations of motion determined from the effective 



Hamiltonian ([34]). This effective Hamiltonian is moti- 
vated from the construction of the quantum Hamiltonian 
operator in a somewhat heuristic fashion. In principle, 
however, there can be additional modifications arising 
from the quantum theory. To properly justify the anal- 
ysis here, more work is required to analyze the quantum 
dynamics by constructing semi-classical states and evolv- 
ing them with the quantum difference equation. The 
accuracy of the effective equations of motion has been 
established in the isotropic case with a massless scalar 
field, and thus we expect that they should give an accu- 
rate representation of the dominant corrections arising in 
the quantum dynamics. 

In addition, we have ignored the inverse triad effects 
arising in loop quantum cosmology. One might ask if 
we do include these effects, does the bounce picture still 
hold. Yet, we can answer this by noting that the inverse 
triad modifications tend to suppress the matter energy 
density. Thus if the suppression is large, the universe be- 
haves more like a vacuum Bianchi I model. Our analysis 
indicates that the bounce still occurs in this case. There- 
fore, the inverse triad effects are not expected to remove 
the bounce, and only change the quantitative behavior. 

We have not considered in this paper additional 
anisotropic Bianchi models. It would be interesting to 
extend the analysis to these models. In particular the 
results in principle could be extended to the Bianchi IX 
model. An interesting question would be whether the 
quantum effects tame the classical chaos in the model 
and again whether a bounce is predicted. Additionally, 
since the Bianchi IX model has been conjectured to be 
of relevance to singularities in general, the analysis may 
provide hints as to what role loop quantum effects play 
with regards to general singularities. 
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APPENDIX A: DETAILS OF THE CLASSICAL 
SOLUTIONS 

In this appendix we discuss in more detail the solutions 
to the classical equations of motion. We show explicitly 
how the universe near the singularity is dominated by 
Kasner like dynamics and how the universe can isotropize 
far away from the singularity for matter with w < \. 

Let us start with the relation (|18p . which, in addition 
to l[T9|). implies 

Pici = K-ih[}Ci + f(t')] (Al) 
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FIG. 1: Mean scale factor a{t) and directional scale factors aj{t) for vacuum pM = case. The initial conditions are chosen 
corresponding to an initially contracting universe. 
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and the time-dependent part: 



fit') = 
1 



(A6) 



±- 



(/Ci +IC2 + /C3) 



2 , 5A{piP2P3) = 



1/2 



We can scale the constants /C/ = ICkj such that (|A5[) 
gives 



FIG. 2: Shear term for vacuum case. Classically, the value is 
a constant, though quantum mechanically it is not constant 
near the bounce. The post-bounce value approaches the pre- 
bounce value, and nowhere does it blow up. 

with the dimensionless constants of motion /C/. Note 
that 

Kh{)Ci - )Cj) = Voaij. (A2) 
We assume that the matter density is in the form 

PM ^ A{piP2P3)-^'+^^/^ (A3) 

with A a constant and w the state parameter. 

The Hamiltonian constraint H = Q with H given by 
(fn)) then yields 

3/(<')' + 2(/Ci-K/C2 + ^3)/(i') (A4) 
+IC2IC3 + /C1/C3 + /C1/C2 = {Kh^y^A {piP2P3)~ , 
which gives the time-independent part:^ 

/C2/C3 + /C1/C3 + /C1/C2 = (A5) 



We can always absorb an arbitrary constant to f{t') and thus 



Kl + K2 + K3 = 1, + '«2 + '^3 = 1' i'^'^) 

which coincide with the "Kasner condition" satisfied by 
the parameters k/ used for the vacuum Bianchi I solu- 
tions (Kasner solutions). [For the Kasner solutions, apart 
from the trivial solution (Minkowski spacetime), two of 
Kj must be positive while the other negative, giving the 
universe expanding in two direction and contracting in 
the other (or the other way around if /C < 0).] Further- 
more, note that the opposite choice of the sign ± in (|A6|) 
amounts to the changes: f{t') —f{t') and JCi — /C/ 
simultaneously, which correspond to the time reversal. 
Therefore, without losing generality, we can stick with 
positive sign for A > 0; ior A < 0, on the other hand, 
± flips sign when the part inside the square bracket of 
(jA6p approaches zero.^ We assume A > in this paper 
for ordinary matter with positive energy. With (jAip and 



change the dichotomy between the time-independent and time- 
dependent parts. However, we pick the particular choice as in 
l lASI I and IIA6I 1 in order to relate /C/ to the standard parameters 
used in the Kasner solutions. 

For A < 0, ± changes signs and thus the solutions of IIA8II may 
encounter recoUapse. This is the case for the model with negative 
cosmological constant A < 0. 
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FIG. 3: Mean scale factor a{t) and directional scale factors ai{t) for a massless scalar field with momentum = lOlp. Initially, 
all directions are contracting and a bounce occurs leading to expansion in all directions. 
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FIG. 4: Shear term for massless scalar field. As in the vacuum 
case, the pre-bounce value is conserved after the bounce. 



H), the equation of motion gives 
1 dpi 1 



Pi dt' 7 



{C2P2 + C3P3) = Kh [/C2 + /C3 + 2/(t') 



(A8) 



(/Ci + /C2 + /C3) 



2 , 3A(piP2P3) ' 



1/2 



For w < 1, the second term in the square bracket of 
is neghgible when the solution approaches the big 
bang singularity; on the opposite side, it becomes domi- 
nant in the large universe limit. ^ Therefore, in the vicin- 



For w = 1, the second term becomes constant and so does f(t')\ 
as a result, the solutions have qualitatively different features. 
This is the case with massless scalar field. In particular, the 
scalar field can be treated as "internal time" and the Kasner 
condition \M} is modified such that "Kasner-unlike" solutions 
(namely, expanding/contracting in all three directions) are also 



ity of the singularity, (|A8|) is approximated as 

1^«^;,(/C2 + /C3), (A9) 
Pi dv 

yielding the solutions very close to Kasner solutions, 
which are highly anisotropic. On the other hand, in the 
large universe limit, (|A8p have the asymptotic behavior 



1 dpi 
pi dt' 



(PlP2P3)^/^ dpi 



Pi dt 
giving the asymptotic solution 



pj{t) cx 



t3(l + ™) 



2^-^ {P1P2P3) 



w ^ -1, 
w — —1. 



(AlO) 
(All) 



As the universe approaches the asymptotic region, the 
three directions are all expanding with the same rate; 
that is, with the inclusion of matter with equation of state 
w < 1, the universe isotropizcs in the expanding phase. 
When the contribution from matter sector is negligible 
and the evolution is essentially the same as the Kasner 
solution as given by (jA9p . we call it "Kasner phase". 
On the opposite, when the matter sector dominates and 
the universe is isotropized as given by (|A10[) . we call it 
"isotropized phase". The situation in between is called 
"transition phase" . 

As anisotropy is concerned, we study the Hubble ra- 
tios, which give the asymptotic behaviors: 



El 

Hj 



di/ai 
cij/aj 



ki/kj 
1 



for a 
for a 



0, 

00, 



(A12) 



which are implied by equations (|A9p and (jAlOp . These 
ratios approach unity in the isotropized phase and ap- 
proach fixed constants in the Kasner phase. ^ 



allowed Il3|l : plus, the anisotropy persists in the expanding phase 
if no other matter content is included. 

Note that, however, in the special case of vacuum (A=0; i.e Kas- 
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FIG. 6: Shear term for radiation. The post-bounce value does 
not equal the pre-bounce values although the value is finite 
through the entire evolution. 



APPENDIX B: DETAILS OF THE EFFECTIVE 
LOOP QUANTUM SOLUTIONS 

In order to affirm the assertions in Section IIIII for the 
generic cases with arbitrary matter, this section deals 
with the detailed analysis for the effective loop quantum 
solutions. The effective dynamics with quantum correc- 
tions is governed by ([551) . (EHl) and which are compli- 
cated to solve analytically for the generic case. However, 
we can still study two extreme cases: a — > oo and a — *■ 0. 

In the limit a — > oo, if /I/C/ — > 0, the Hamilton's equa- 
tions plus the Hamiltonian constraint simply reduce to 
their classical counterparts and the effective solutions are 
virtually the same as the classical one. We do not know 
a priori whether fljci —^ for large a and whether the 
quantum effect does not spoil the semi-classicality. But if 
the classical solutions without quantum corrections leads 



ner model) or of scalar matter {w = 1), the Hubble ratios Hj/Hj 
are constant throughout the entire evolution. The classical solu- 
tion does not isotropize and remains in the Kasner phase. 
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FIG. 7: Ratios of directional Hubble rates for radiation. At 
late times, the ratios approach unitary indicative of equal 
expansion rates in all directions as the universe isotropizes. 



to small value of fijcj for large universe, we can con- 
clude that the effective dynamics admits the solutions 
with semi-classicality for large a. To check this, consider 
the classical solutions, which have 

PICI = K-,h[)Cl+fit')]^K-fhf{t') 



» Ty ipiP2P3)~ 

for a ^ CX3 by (jAT|) . We then have, by (|Allj) . 



ci (X a 



fj^icj cx a 



3(1 + ™) 



(Bl) 



(B2) 
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This shows that the quantity /2/C/ of classical soultions 
is decreasing to zero in the large universe limit for w > 
— 1. Therefore, for — 1 < < 1, at large universe, the 
semi-classicality can be retained and the loop quantum 
corrections are indeed negligible for those solutions which 
are semi-classical at late times. The evolution simply 
follows the classical trajectory for the large universe.^ 

Now, let us study the other extreme as the universe 
approaches the classical singularity (a — > 0). In the 
backward evolution, before the singularity is reached, 
the quantum corrections will take effect and the big 
bounce is expected to take place. Therefore, at some 
point, cos(/i/C/) vanishes and flips sign. Assuming that 
cos(/i/C/) in different directions flip sign at only slightly 
different moments, we then have cos(/2/c/) ~ cos{p,jcj) 
0) in the vicinity of the big bounce. With this ap- 
proximation, close to the epoch of the big bounce, pS]) 
yields^ 



(B3) 



1 — w 



K7 [cos(^/c/) - cos(/ijcj)] (P1P2P3)' 



It then follows from (El) that 



sin(^/C/) « Kjh' 



Ai/2 

3/2 
Pi 

3/2 ■ 
Pi 



ICi 



A 



{P1P2P3) " 



provided 



ICi > 



A 



(P1P2P3)" 



(B6) 



(B7) 



Taking (jB4p into (|55|) and expressing cos a; = ±(1 
sin^a;)^/'^ with the help of (|B6p . we get for instance 



— « ±kH{IC2 + IC3 
Pi at' 



Pi 



{JCi + 



1/2 



±Kh {K.2 + /C3) 



Pl.crit 
Pi 



1/2 



(B8) 



Thereby, we can write 

sin(^/c/) 
Qi =Pi : 



(B4) 



in accordance with the classical counterpart (|Aip ; conse- 
quently the Hamiltonian constraint (|39p gives the same 
f{t') as given by (jA6p . To proceed further, in the fol- 
lowing, we consider three cases separately: quantum cor- 
rections take effect (i) in the Kasner phase; (ii) in the 
isotropized phase; and (iii) in the transition phase. 

In Case (i), the first term in the square bracket domi- 
nates over the second term in (jA6P : thereby. 



where the critical value pi^crit is given by the Planck 
length square £p^ times a numerical factor: 

« 19.0l7|/C/|^/^4i- (B9) 

Therefore, the big bang singularity is replaced by the 
bounces whenever each of pi approaches its critical value 
P/.crit- The bounces occur up to three times, once in each 
diagonal direction. If we define the directional density gj 
in the /-direction as: 



Qi 



ip] ' 



(BIO) 



f{t') 



IC 
3" 

A 



JC 
3" 

{P1P2P3) 



sA{piP2P3y 



(B5) 



the above statement can be rephrased to say: The big 
bounces take place whenever each of the directional den- 
sities reaches the critical value 



Pcrit = 3(k7^A)" 



.82ppi. 



(Bll) 



* This is also true for the special cases of vacuum and of scalar 

matter with ui = 1, in both of which picj = constant and thus 

_ —3/2 

/i/C/ oc pj — > classically. For the case of the cosmological 
constant {w = —1), the cosmological constant A has to be very 
small to admit the semi-classicality. 
^ The assumption that cos{/i/C/) flip signs at only slightly differ- 
ent moments could be wrong in general. However, note that 
even if cos(p,jCj) ^ cos(/2jCj) does not hold very well, the van- 
ishing of IIB3I 1 is still a good approximation near the bounce, 
since {pip2P3)^^~^"' ^"^ = a^^^~^"' — » for ui < 1 when a is small 
enough. (But the condition cos(/i/c/) cos(/2jcj) makes the 
approximation | |B3| I even more accurate.) In the special case for 
w = 1, l|B3| l vanishes exactly simply because the factor (1 — w) /2 



Note that gi have the same dimension as pM and more- 
over we have 



pIqi +pIq2 +Ph3 



with 



E2 



18 



("12 + " 



a 



317 



(B12) 



(B13) 



being identical to the shear factor of the classical solu- 
tion. This suggests that gj can be roughly interpreted as 
the "energy density carried from the classical anisotropic 
shear portioned for the /-direction" (and thus the name) . 
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To meet the condition (|B7p of Case (i), we take the 
critical values p/.crit hito (|B7p and find the criterion for 
Case (i) to be 



IK1K2K3 



(B14) 



Next, in Case (ii), on the other hand, the second term 
in the square bracket of (|A6p dominates and we have 



fit') 



A 



It then follows from (jBij) that 
Ai/2 



sin(/x/c/) 
provided 



A 



p3/2 V 3Kft2 



(P1P2P3) 



With the help of (|B6| again, ^ gives 
1 dpi 



pi dt' 



±2Kh 



A 
3^ 



(B15) 



(B16) 



(B17) 



(B18) 
1/2 



(P1P2P3)" 



The big bang singularity is again replaced by the 
bounces and the bouncing points of pi are roughly equal 
in all three direction; i.e., each pi is bounced when 
Pi ~ Pcrit with 



Pcrit : = 



(B19) 



For more generic cases (with multiple matters), this 
means that the bounces take place near the point when 
the (total) matter density approaches its critical value 



Pent — ^Pcrit 



(B20) 



begin the same as the critical value given in (IBlip . 

To meet the condition (jBlVp of Case (ii), we take the 
critical value pcrit into (|B17p and find the criterion for 
Case (ii) to be 



A < 7"'-^/C) 



(B21) 



which is exactly the opposite of the criterion (jB14p if we 
have Ki ~ 0(1) • 

Finally, in Case (iii), we have 



which yields 



P/,crit ~ Pcrit- 



(B22) 



(B23) 



Therefore, in Case (iii), the big bang singularity is re- 
placed by big bounces as well and the bouncing points 
of Pi are between p/,crit given by (|B9p and pcrit given by 
(jBigp (p/,ciit and pcrit are now in the same order). 

To sum up, the big bang singularity is replace by big 
bounces due to the fact that the gravity with loop quan- 
tum corrections becomes repulsive at some point when 
the universe is near the singularity. In Cases (i), this 
happens when the directional density gi reaches Planck- 
ian energy density .86ppi; in Case (ii), it happens when 
the matter density pM approaches ~ .86ppi. In all cases, 
the big bounce takes place whenever either of gi or pM 
approaches 0{pp\) first. In a sense, there is competi- 
tion between the energy density carried from the classical 
anisotropic shear (gi) and the matter density (pm) to be 
the indication of occurrence of the big bounces. 

Also notice that as mentioned in Section lllli 
(post bounce) 7^ (pre bounce) in general due to the 
fact that the anisotropies are smeared by the quantum 
corrections. However, in Case (i), the information of 
anisotropies persists during the bouncing period and we 
shall have S^(post bounce) « E^(pre bounce). This can 
be seen from (IBSp . in which the evolutions for different 
Pi are decoupled around the bounces and therefore the 
constants /C/ (which dictate anisotropies) are unchanged 
before and after the bounce. 

The equations of motion given by l|35p and (|36p can be 
solved numerically once the initial values pi{t = to) and 
ci(t = to) are given. Equivalently, for given pi{to), spec- 
ifying JCi is an alternative way to specify c/(to). [Given 
f{to) can be obtained via (|A6p and then c/(io) are 
fixed by (|Aip if the initial condition is in the classical 
regime and by (jB4p if in quantum regime.] By changing 
the parameters A and JC, we can get the three different 
cases as discussed above. The numerical solutions in the 
presence of radiation are solved in terms of proper time t 
and depicted in Figures [8| [9l and [TOl for these three cases 
respectively. 



APPENDIX C: ALTERNATE QUANTIZATION tion|TTTl The major difference is that the behavior of the 

parameters pi varies with the triad components differ- 
In this section we consider the effective dynamics of ently Specifically we take the dependence as p/ cx 1/a/ 
an alternative quantization scheme mentioned in Sec- 
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FIG. 8: Case (i): Quantum corrections take effect in the Kasner phase. With w = 1/3 (radiation filled); ki — —2/7, 
K2 = 3/7, K3 = 6/7; /C = 1. X 10^ A = 1. x lO^fi^T"^ and pi(to) = 3. x 10*41, P2(to) = 2. x 10*4i, P3(to) = 1. x lO'*^!- 
(Also set 7 = 1.) The red curves are for pi, ai; green for p2, a.2', and blue for ps, 03. Solid lines are the solution to the effective 
loop quantum evolution; dashed lines to the classical evolution. The proper time t is offset such that the classical singularity is 
at the origin of t. The values of pi,crit are pointed by the colored arrows and Pcrit by the black one. In this case, the bouncing 
point of each pi matches p/,crit very precisely and we have p/,crit S> Pcrit (or say pM Qi ^ ppi near the bounces). [Note that 
the bounce of pa is out of the shown range. The isotropized phases on both sides of the bounces are also out of plot.] 
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FIG. 9; Case (ii): Quantum corrections take effect in the isotropized phase. With w = 1/3 (radiation filled); 
Ki = -2/7, K2 = 3/7, = 6/7; /C = 1. x 10^ A = 1. x lO^fi^T"^ and pi(to) = 9. x W^i, P2(to) = 6. x 10*4i, 
P3(io) = 3. x 10*4i- In Ihis case, the bouncing point of each pi roughly matches pcrit and we have p/,crit ^ Pcrit (or say 
Qi ^ Pm ~ Ppi near the bounces). [Note that in the backward evolution the contracting curve of the classical ai eventually 
becomes expanding at the epoch extremely close to the singularity, indicating that the Kasner phase does occur classically 
although almost invisible in the plot.] 
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FIG. 10: Case (iii): Quantum corrections take effect in the transition phase. With w — 1/3 (radiation filled); 
Ki = -2/7, K2 = 3/7, K3 = 6/7; /C = 1. x 10^ A = 2. x 10''ft4i""^ and pi(to) = 3. x 10''4i, P2{to) = 2. x 10*4i, 
P3{to) = 1. X 10*4i- In this case, we have p7,crit ~ Pcrit, roughly around which all pi are bounced (or say gi ~ pM ~ ppi near 
the bounces). 



14 



with the exact behavior for /ii given as 



Ml 



P2P3 



(CI) 



The effective Hamiltonian is similar in form to that in 
([34|) and is given exphcitly by 

1 f sin(^2C2) sinifi^cs) 
K7 I M2M3 

+cychc tcrmsj + P1P2P3 Pm, (C2) 

which becomes 

rtcff = TT-"^ Sm(^2C2)sm(^3C3) 

+cychc tcrmsj + P1P2P3 Pm- (C3) 

With this, the vanishing of the constraint leads to the 
simple relation 



1 



Pm 



sin(/i2C2) sin(/i3C3) + cyclic terms 



and from this wc can deduce an important observation. 
Namely, since all of the sin terms are bounded we find 
that the energy density is also bounded as 



Pm < 



Merit J 



(C4) 



where pcrit is numerically the same value as in the 
isotropic case 



Merit 



.82p 



PI- 



(C5) 



The fact that pM is bounded implies that the big-bang 
singularity is resolved since classically the energy density 
blows up there. Note that any bounce does not neces- 
sarily occur when pM = Merit and we shall see, that with 
anisotropics the bounce occurs at lower energy densities. 

We can derive Hamilton's equations in the same fash- 
ion as in Section IIIII and it can be shown that 



Pici - pjcj = 7Vb aij, 



(C6) 



where ajj are constants of motion. Note that this ex- 
act relation was satisfied in the classical case in ([T9|. 
However, with this effective Hamiltonian, c/ are not as 
simply related to dj as in the classical case given in ([T5|) . 
The full relation is to be determined from the Hamilton's 
equation for pi which give for instance 



Pi 



Pi 



A-f 



cos(/ziCi) sin(^2C2)+sin(^3C3) . (C7) 



Because of this, the shear term = jg[{Hi ~ H^)^ -t- 
(H2 — -ffs)^ -t- (-ffs — i/i )^] is no longer constant. However, 



in the classical limit where /l/c/ ^ 1, we can calculate 
the shear term to be Y? k, j^{ci\2 + 0^23 + '^si) in the 
classical case. 

Therefore if we consider the behavior of the shear term 
and begin with a semi-classical collapsing universe, we 
find that initially Y? is constant with value 



18 



12 



*23 



^3U 



(C8) 



Through the bounce S is not constant, but after the 
bounce in the expanding phase where once again /2/c/ ^ 
1, we find that the shear takes on the same value. This 
is because equation (|C6p holds throughout the evolution. 
Therefore it is a rather general feature of this effective 
Hamiltonian, that the shear term E is conserved before 
and after the collision. Note that this is in contrast with 
the results presented in the body of this paper, where E 
was only conserved for a massless scalar field. However, 
both quantizations imply that the shear term is bounded 
and the anisotropics do not blow up. 

Because of the complexity of the equations of motion, 
it is highly non-trivial to derive a generalized Freidmann 
equation. We can however if we expand the effective 
Hamiltonian to second order in /2/C/. Then equations 
(|C6P and (|C7|) can be used to derive an approximate gen- 
eralized Friedmann equation given by 



a = oMM 1 

J ^ Merit 



E2 
^6 



a^pcrit Kal2p^j.ij 



+ 0((/2,c/)4). (C9) 



If we solve this for the matter energy density at the 
bounce [H = 0) we find 



Mbo 



1/2 Perit g 



(Merit 



) (Merit 




, (CIO) 



which implies that the energy density is bounded below 
Pcrit in accordance with the prediction from the vanishing 
of the effective Hamiltonian. Note that is a constant 
of motion and only related to the shear term E^ in the 
classical limit p,jcj <C 1. Note also, that this generalized 
Friedmann equation is valid in the nearly isotropic limit 
which can be understood from the relations (|C6P which 
imply 



M/c/ - Pjcj 



jVAaij 



(Cll) 



If the right hand side is not small (i.e., when the aniso- 
topies are large), then at least one /2/C/ is guaranteed to 
be large and hence the generalized Friedmann equation 
would need to be calculated to higher order to provide a 
good approximation. 
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We note that both the quantization scheme presented 
in this appendix and that given in the body of this work 
agree as to the quahtative nature of the bouncing uni- 
verse. Namely, a bounce still occurs with the inclusion 



of anisotropics and that the shear term does not blow up 
implying that the anisotropies remain finite through the 
bounce. 
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